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The layered perovskite compound NaV2 05 exhibits charge 
ordering (CO) and spin gap formation which is driven by 
inter- site Coulomb forces and exchange dimerisation. Ex- 
perimentally it was found that the CO transition is rapidly 
reduced by doping. Charge ordering in NaV^Os has a quasi 
ID character due to inter- ladder frustration in the Trellis lat- 
tice. Based on an extended pseudospin (T=l) approach for 
hole doped NaV^Os ladders we construct an Ising like model 
that explains qualitatively the destruction of the CO phase 
as function of hole doping S and show how the CO parameter 
depends on 5 and the ratio of Coulomb and hopping energies. 

PACS numbers: 71.I0.-w, 71.27.+a, 71.28.+d 



I. INTRODUCTION 

The layered perovskite NaV20s shows an interesting 
interplay between charge ordering (CO) of 3d- electrons 
on the rungs of V-V ladders and a spin- Peierls tran- 
sition associated with spin gap formation first observed 
in 1 below T c = 34 K. This transition is associated with 
a distortion of the high temperature orthorhombic struc- 
ture of NaV205 . The latter can be viewed as Trellis 
lattice with ab- planes stacked along c. They consist of 
ladders with V-O-V rungs (inset of Fig. 1) shifted alter- 
natingly along b by half a lattice constant. The precise 
low temperature structure is still controversial 2-6 with 
the proposed models differing in the number of inequiva- 
lent V and Na sites and an associated charge ordering on 
each or only on every second ladder. The slow progress 
and ambiguity of experimental results has hindered the- 
oretical understanding considerably. However it has be- 
come clear now that there is a close connection between 
charge ordering on the V-V rungs and exchange dimeri- 
sation on the ladders 7 which leads to the spin gap for- 
mation directly seen in inelastic neutron scattering 8 , ft 
was realized early on 9 that the CO on the V- rungs may 
be described within an Ising model in a transverse field 
(ITP) where the isospin (S=|) represents the charge de- 
gree of freedom on the V-V rung with inter- site Coulomb 
forces correspond to Ising interactions and the intra-rung 
hopping to a transverse field. This is due to the charge 
transfer (CT) gap of - 0.4-0.5 eV 11 ' 12 which leads to an 
insulating state of NaV20s despite having only one 3d- 
electron per rung (quarter-filled ladders) . Therefore only 
singly occupied rungs have to be considered which may 
be represented by an isospin and the real spin. A detailed 
model of the combined Ising- spin Peierls system and its 
phase transitions has been discussed in 7 . 



The charge ordering in a Trellis lattice compound is 
severely hindered by frustration of inter-site Coulomb 
forces which is directly obvious in the Ising model rep- 
resentation. In fact on a rigid lattice CO would be pro- 
hibited since the lattice of V-V rungs corresponds to a 
trigonal covering lattice of Ising spins. The simultane- 
ous structural distortion at T c reduces the frustration 
by leading to inequivalent adjacent ladders. The driving 
mechanism for the transition are inter- site Coulomb in- 
teractions and exchange-striction. The gain in Coulomb 
energy (Vb) by complete CO in a rung is in competition 
with the loss of intra-rung kinetic energy (t a ), the control 
parameter for CO is therefore A =2t a /Vb (for convenience 
the inverse definition of 13 is used here). The very low 
CO transition temperature suggest that in NaV20s one 
is close to the quantum critical point A c =l of the cor- 
responding ITP model 13 . The charge ordering is then 
essentially of ID Ising type, this means one has long 
range correlations along the ladder and no correlations 
between them, this was indeed obtained in exact diago- 
nalisation studies of the ITP 14 . The ordering may then 
be viewed as arising from Ising like correlations along the 
chain which lead to a true 2D broken symmetry state at 
finite T due to a staggered longitudinal Ising field set up 
by the distortion of the neighboring ladders 13,7 . 

In this context it is suggestive that even small per- 
turbations of the ID Ising correlations might reduce or 
destroy the CO state in NaV20s . This can be achieved 
by reducing the filling factor n of the ID ladders be- 
low j by doping with holes. This introduces "empty" 
rungs into the ladder which cut the Ising bonds. The 
ensuing destruction of long range ID correlations should 
then strongly reduce T c as function of the hole concentra- 
tion S (n=j-S). This has indeed been found by introduc- 
ing holes in the ladders through Na- deficiency doping 18 
where a few per cent holes are sufficient to destroy the CO 
state and the associated spin Peierls transition. Rapid 
suppression of charge order has also been found in var- 
ious other doping series, i.e. replacing Na by Li and K 
(isoelectronic) or Ca (electron-doping). 15-17 This may be 
observed directly by specific heat measurements 16 which 
show a progressive suppression of AC(T C ) with increas- 
ing doping. It is also seen in the susceptibility 18 ' 15 which 
shows a closing of the CO induced spin gap. 

Most importantly in this doping range NaV20s re- 
mains an insulator which is not easy to understand within 
a Hubbard like model for the quarter filled ladder 5 ' 7 , pos- 
sibly ID loalization effects play a role, indeed the con- 
ductivity was found to exhibit variable range hopping 
behaviour 18 for hole doping. The notion that charge or- 
dering T c is reduced because doping creates empty rungs 
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which cut the Ising bonds along the ladder is however 
completely static and invokes a disordered state. It lacks 
the important ingredient of kinetic energy gain due to 
the hopping of holes between adjacent rungs leading to 
a homogeneous state. This problem has sofar not been 
considered theoretically. 

In this work we will give a more solid theoretical basis 
for the influence of hole doping on CO in NaV205 . This 
will again be done in the context of an effective isospin 
approach. In sect. II we give the basics of the electronic 
structure of NaV 2 5 and discuss the isospin model for 
the undoped case in sect. III. Subsequently in sect. IV we 
show that even in the doped case the low energy Hamil- 
tonian may be projected to an effective Ising type form, 
albeit with an Ising spin T=l to accomodate the holes 
on empty rungs. In sect. V we discuss how order param- 
eter, kinetic energy and chemical potential (doping level) 
are represented in the isospin formalism and their respec- 
tive coupled mean field equations in the ground state in 
sect. VI. They are solved in sect. VII numerically and we 
discuss the variation of the CO parameter with the con- 
trol parameter A and the hole doping S which leads to the 
zero temperature CO phase diagram. Finally sect. VIII 
gives summary and outlook. 



II. ELECTRONIC STRUCTURE OF NaV 2 O s 

The electronic structure of the charge transfer insula- 
tor NaV205 has been investigated 10 ' 11 within the LDA 
and LDA+U approach respectively. A mapping of the re- 
sulting band structure to an extended tight binding (TB) 
model involving both V 3d and O 2p orbitals has also 
been performed 11 and an excellent fitting of the LDA+U 
bandstructure was obtained. For our present purpose it 
is sufficient to use a reduced TB model involving only 
V 3d xy orbitals. The hopping elements involved in this 
model are defined in Fig. 1. There are four V atoms 
in the unit cell resulting in two 3d xy occupied valence 
and two empty 3d xy conduction bands separated by a 
charge transfer gap of about 0.3 meV which is largely 
determined by the effective intra- rung hopping t a . In 
the effective TB model this gap may be interpreted as 
bonding- antibonding gap of V 3d xy orbitals of a rung. 
Furthermore the valence and conduction band are each 
almost twofold degenerate due to the weak hopping be- 
tween ladders. The 3d xy bonding bands have a consid- 
erable dispersion of order tf,+t<j along the b- direction 
while the dispersion t^-t^ of 3d xy antibonding bands is 
almost zero due to accidental near equality t;, ~ td- 



III. THE ISOSPIN MODEL HAMILTONIAN FOR 
UNDOPED NaV 2 5 

The basic many body Hamiltonian de- 
scribing NaV 2 05 is of the single band {3d xy ) extended 



Hubbard type 9 ' 7 . As shown in 9 , due to the presence of 
a CT gap the charge degrees of freedom may be repre- 
sented by an effective ITP Hamiltonian (T=^) which is 
obtained by projecting the original Hamiltonian to the 
subspace of singly occupied rungs 7 . The states | — 5), 
1 1) of the T=| isospin describes occupied d xy - orbitals 
on the left and right V atom of each rung respectively. A 
simplified version of the resulting ITP Hamiltonian for a 
single ladder is given by 13 

Hitp — 2V;, y^T^Tzj + -] - 2t a Tix 

(ij) « 

- h ° s J2(-iyT zi (1) 

i 

Here t a is the hopping across the rung and Vf, the 
Coulomb interaction along the leg of a ladder. The last 
term is due to a staggered longitudinal isospin field h s j 
~(-iyh^ which simulates the influence of the neighbor- 
ing ladders 13 . For A = 2t a /V& < A c =l one obtains an 
AF ordering of isospins which corresponds to the zig-zag 
ordering of 3d- electrons on the ladder as illustrated in 
Fig. 1. It was shown in 7 that the interaction parameters 
in Eq. I are renormalized by virtual processes via dou- 
bly occupied rungs or sites which depends on the spin 
state. This leads to a coupling of isospin (charge) and 
spin degrees of freedom whose contribution is neglected 
here. 

The zig-zag CO in a ladder is characterized by the 
order parameter (T Z i) = (— l) l a; at zero temperature with 
(A<1) 



-1/3 



(2) 



The exact solution of the ID ITF model has the ex- 
ponent /? = I and A c =l. It is understood that an in- 
finitesimal staggered field h S i has to be applied to lift the 
twofold degeneracy of the ordered ground state. A trun- 
cated mean field (mf) approximation suppressing double 
counting of bonds in the molecular field leads to a criti- 
cal value A c =l. This is the same quantum critical point 
that separates the CO (A <f ) phase from the disordered 
phase (A >f) as for the exact solution 13 . Naturally the 
mf critical OP exponent /3 = | is different. 



IV. T=l ISOSPIN MODEL FOR DOPED NaV 2 5 

The charge ordering described by the ITP model and 
its implications for the (spin-) superexchange have been 
described in detail in 13 ' 7 . In the present work we want 
to investigate whether the convenient and physically ap- 
pealing isospin approach for CO in NaV2 0s can also be 
applied to the hole doped compound to study the sup- 
pression of CO. For this purpose the hole states, that is 
empty V-V rungs with zero 3d- occupation, also have to 
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be represented by an isospin degree of freedom. As in 
the stoichiometric case 13 ' 7 doubly occupied rung states 
either have a large charge transfer energy (t a ) or a large 
on-site Coulomb repulsion and therefore contribute only 
as intermediate states by renormalizing the interaction 
parameters of the effective Hamiltonian. Thus the model 
has to incorporate three basis states, two singly occu- 
pied and one empty rung state which suggests to use a 
T=l isospin instead of T=^ for the stoichiometric (<5=0) 
case. The mapping of electronic states to isospin states 
for rung i is then defined as 



-1), 
0), 



|10>i 
|00)i 



4l00)i 



(3) 



\T Z = 1), = |01)i = 4|00), 



where P +Pi = l and Pi=P L +P fl holds. The other 
projector properties PqPi= PiPo=0 and P^= Pa 
(A=0,L,R) are easily verified. Using these projection op- 
erators the intra- ladder Coulomb interactions V& be- 
tween adjacent rungs i,j may be represented by the 
Hamiltonian 



H c (i,j) = \(V b 



V d )T 2 zi T 2 zj 



Vd)T zi T Z j 



(5) 



Compared to the T=| case of Eq. (1) He also con- 
tains now biquadratic terms in the isospin. Similarly we 
construct the intra- rung kinetic energy responsible for 
the charge transfer gap. Since the R,L- occupation of 
each rung is now described by \T Z = ±1) respectively the 
intra-rung hopping term must change the isospin by ±2 
units. The corresponding rung Hamiltonian then reads 



> 

CD 




FIG. 1. Mean field eigenvalues (full lines) E K (k= 0-2) of 
Eq. (14) and 'chemical potential' fi (dashed line) as function 
of the control parameter A. Arrow indicates CO transition. 
Inset shows the NaV20s ladder along b with effective V-V 
rungs along a. Closed and open circles correspond to occupied 
and empty d xy - orbitals. Hole doped rung has two empty 
d xy - orbitals. Only the ground state (fc=0) is occupied for 
each rung with a hole count given by Eq. 17. Hopping and 
interaction parameters are indicated. We used t a = 0.380 eV, 
tb= td =0.170 eV; for numerical reasons the latter is twice 
the value used in 13 . 

with |00)j denoting the empty rung state. We first ex- 
press the inter- rung Coulomb energy already present in 
Eq. (1) in the above extended isospin basis. The projec- 
tors to this basis may be written in terms of the T=l 
isospin operators T= (T^ 



,Tj,,T z ) as 



1, 



Pl = ^T Z {T Z + 1) 

Pr = \t z {t z 1) 

Pn = 1 - T? 



(4) 



H K {i) 



\t a {T 2 +l + T 2 _ 



MT 2 X 



rp2 



(6) 



which should be compared to the second term in 
Eq. (1) for T=i. Finally there has to be an additional 
term describing the hole kinetic energy, i.e. when the 
empty rung (Fig.l) is shifted along the ladder direction. 
Such processes are described by inter- rung terms like 

H' K (i,j) = 2t + [O zy {i)0\ y {j) + Ot y (i)O zy (j)] 

+ 2t-[O zx (i)0\ x {j) + Ol x (i)O zx {j)\ (7) 



with i,j denoting rungs along the ladder. Here 
t±=2 (tft±td) and Z j / =iT 2 Tj ) and O zx = r T z r £ x may be 
interpreted as a quadrupolar operator in the T=l isospin 
space. In the following we consider only the special case 
tb ~ td which is appropriate for NaX^Os 11 . Then the 
second term above vanishes and the total T=l model 
Hamiltonian for the doped ladder now reads 



H - ^2\ eT zi + ta(T xi -Tyi)- h° si T z , 

i 

+ 2t b Y,[O z yOH+0%Oiy] 



(8) 



(H) 



Here e is the on-site orbital energy and h°(i)= hg(— 1)' 
is again a staggered field due to the influence of neigh- 
boring ladders, (ij) denotes n.n. rungs along the ladder. 
Furthermore we have defined V± = \{V b ± V d ). If we ne- 
glect the ladder diagonal Coulomb interaction (Vd =0) 
then V± = ^Vfc leading to a single control parameter A= 
2t a /V&. Although we expect the isospin model a conve- 
nient framework to discuss the energetics of charge order- 
ing as function of doping it can however not give answer 
to the question why the doped system stays an insulator, 
for this problem one also has to include the doubly oc- 
cupied rung states. It was argued in 5 ' 7 that even under 
a finite charge transfer between ladders, i.e. deviations 
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from quarter-filling, an insulating state can be expected 
for the relevant hopping and interaction parameters re- 
alised in NaV2C>5 . 



A x 0.5 
h- 





5=0.0 
0.03 
0.06 
0.09 


0.09 
0.06 

8=0.03 


1 



0.2 



0.4 



0.6 



0.8 



*=2L/V h 



FIG. 2. zig-zag charge order parameter (T z ) as function 
of control parameter A for small and moderate hole doping 
S (full lines). Dashed lines represent average hole hopping 
energy (T x ) (units: t 6 ), for <5=0: (T x ) = 



V. CHARGE ORDER PARAMETER, DOPING 
AND 'CHEMICAL POTENTIAL' 

It is obvious that the model defined by eq. 8 is much 
more complex due to the additional hole state (T 2 =0) 
and has no exact ID solution like the T=i ITF model. 
As mentioned previously the truncated mf approximation 
for the charge order parameter of the ITF model has the 
same quantum critical point A c =l as the exact solution. 
It is therefore plausible that a mf treatment of Eq. 8 in the 
doped case still gives a qualitatively correct picture for 
the conditions of charge ordering. As in the ITF case the 
order parameter for Vb >Vd should be staggered, i.e zig- 
zag CO should be realized, then for the two sublattices 
we have 



(T zA ) = -(T zB ) 



(9) 



Now x not only depends on the control parameter 
A=2t a /V& but also on the number of holes or unoccupied 
rungs. The representation of the hole number operator 
in fermionic and isospin language is given by 



Ai C Li 



in 



CRi 



5i = l-hi= , _ 

5i = l-P iL -P iR = l- Tl (10) 
so that the hole concentration is obtained as 



5 = 1-(T?) 



(11) 



Since it is externally fixed by the level of doping the 
local orbital energy has to be adjusted selfconsistcntly 



for each parameter set such that S is obtained. This can 
be done by adding a 'chemical potential' , i.e. replacing 
e — > e — /i. For each microscopic parameter set and order 
parameter x, the chemical potential fi has to be adjusted 
appropriately to obtain the physical hole number. 




FIG. 3. upper panel: melting of CO with hole doping S for 
various A given in decreasing order. The absolute slope value 
-(dx/d<5) increases strongly when approaching the quantum 
critical point A c =l of CO; this is shown in the inset, lower 
panel: Hole kinetic energy (T x ) (in units of tb) as function of 8 
in the complete doping range, plotted for five values of A cor- 
responding to curves in increasing order. Downward arrows 
indicate CO transition (x=0) where the curves for different A 
merge. 



VI. MEAN FIELD APPROXIMATION AND 
SOLUTION FOR CO 

The usefulness of a simple mf solution relics on the 
Ising type nature of the present model. In this approxi- 
mation the doping term <~ % in Eq. 8 simplifies consider- 
ably using O zy +0\ = T x 



and O zy -0\ y -- 



■iTy, we obtain 
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H' Kmf (i) 



= -2t b [O zy (O yz )+O yz (O zy )} 
= -t b ((T x )T x (i) + (T v )T y (i)) 
= -t b (T) (cos 4>T X + sin0T y ) 



(12) 



where (T)=y/ (T x ) 2 + (T y ) 2 and <j> which is given by 
tan -1 cj)=(Ty)/(T x ) is an arbitrary angle which may be 
chosen as 0=0. Then only (T x ) appears and the quan- 
tity W= t b (T x ) may be interpreted as the average inter- 
rung hopping amplitude of the hole. Together with the 
remaining terms in Eq. 8 we get the following rung mf 
Hamiltonian for A,B sublattices: 



rrmf 



H 



mf 



[e + V + {l-5)-n]T 2 zA -t a {T 2 xA - 

-(h° s +V-x)T zA -WT xA 

[e + V + {l-S)-^]T 2 B -t a {T 2 xB 

+ (h° + V-x)T zB - WT xB 



T yA) 



Tin) (13) 



The effective mf orbital energies and staggered fields 
are then given by 

e = e- n + V+il-S) 
h s = h° s +V-x 

In the T=l isospin space the mf Hamiltonian for A,B 
sublattices can explicitly be written as 



H mf - 

A,B — 




-W/y/2 


-W/V2 




(14) 



where A,B corresponds to upper and lower signs re- 
spectively. These equations have to be selfconsistently 
solved for the charge order parameter x, keeping the 
hole concentration fixed by varying the chemical poten- 
tial [i. Only the T=0 case is considered, then the expec- 



tation values (X) = (tp™ \X\tp™ ) are obtained from the 
ground state i/j™* solution of Eq. 14. For the actual cal- 
culation, in addition to assuming t a ~ td which has been 
justified in , we take for simplicity Vd — 0. If we sup- 
press the hole hopping described by the matrix elements 
~W for the moment, the T z =0 and T z = ±1 subspaces 
disconnect and the mf ground state and its energy of the 
isolated rungs can be written explicitly as 13 

1 



£o=e--V(Ae) 2 + (2t„) 2 
W / )i=«ii|l)+«2i|-l) 



(15) 



where Ae=V;,x is the splitting of L,R orbital energies 
and the staggered mf coefficients are given by 



\ 



1± 



(-IJ'Aei 



v/(Ae) 2 + (2t a y 



(16) 



When the effect of hole hopping along the legs (W> 0) 
is now included the mf eigenstates will be mixtures of the 



isolated rung states an d the hole state |0)j 

hole concentration is then given by 



u m f\ 



The 



(17) 



The new mf- ground state energy, its wave function 
and associatied expectation values of T x , T z and T 2 Z then 
have to be calculated selfconsistently from Eq. 14. The 
numerical results will be discussed in the next section. 



VII. NUMERICAL SOLUTION, DOPING 
DEPENDENCE OF CO AND PHASE DIAGRAM 

In this section the numerical results for the above 
model will be presented. The main goal is the determi- 
nation of the charge order parameter x(S, A) as function 
of the hole concentration 6 and the control parameter A= 
2t a /Vb. To achieve a given hole number 5= 1-(T 2 } the 
'chemical potential' ji has to be adjusted selfconsistently 
together with the determination of the order parameter 
x=(T z ) and the normalized kinetic energy (T x }=W/t b . 
Physically /i corresponds to an orbital energy shift of the 
empty and singly occupied rung states. This energy shift 
together with the mf eigenvalues of Eq. 14 as function of 
the control parameter A is shown in Fig. 1 for (5=0.09. 
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FIG. 4. CO/MV phase diagram in the S, A- plane. Contour 
lines for CO phase with 8={T Z ) for x=0.01, 0.25, 0.5, 0.75 
(moving to lower left corner) are also shown. MV denotes the 
disordered homogeneus mixed valent phase. 

The corresponding variation of the charge order pa- 
rameter x for several hole concentrations is shown in 
Fig. 2. As the hole concentration increases our results 
show the expected qualitative behaviour: i) the CO even 
for large Y b (small A) is no longer complete due to the 
presence of empty rung states and its variation with A 
becomes essentially flat, ii) the critical value A c for the 
CO transition to appear is reduced from the value A c =l 
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of the undoped ITF model. In addition the kinetic en- 
ergy (T x ) is shown. It increases almost linearly with A= 
2t a /V{, and rapidly with the hole concentration. 

The complementary reduction of charge order as func- 
tion of the hole concentration for constant A is shown in 
Fig. 3. Most significantly for small d one observes a rapid 
reduction of the CO- parameter with increasing hole con- 
centration. The slope of this reduction increases even 
further when A approaches the quantum critical point 
A c =l of CO from below. This is shown in more detail in 
the inset of Fig. 3. The origin of this drastic reduction of 
CO for small 6 lies in the rapid increase of hole kinetic 
energy with doping as can bee seen in the lower panel of 
Fig. 3. 

Our results shown in Fig. 3 describe qualitatively cor- 
rect the experimental situation, namely a destruction of 
the CO state in NaV2 05 by a few percent hole doping 16 . 
For example assuming a value A =0.95 close to A c =1 one 
obtains a critical concentration S — 0.03. The maximum 
size of CO at 6=0 for this A- value is already consider- 
ably reduced to x ~ 0.3. The experimental value of the 
nominal hole concentration for suppression of CO in Na- 
deficient NaV 2 5 is indeed about S = 0.03 15 . 

Finally we summarize the results of our calculations in 
Fig. 4 which shows the phase diagram of CO in the (5,A)- 
plane and the size of the charge order parameter x as a 
contour plot in the ordered regime. 

VIII. SUMMARY AND OULOOK 

In this work we have proposed an extended isospin 
approach proposed in 9 to the charge ordering transition 
in Na- deficient (hole-doped) NaV20s . The additional 
hole (empty rung) state can be acommodated within a 
T=l isospin as the T z =0 component wheras the (left 
and right) occupied rung states are represented by the 
T z = ±1 components. The effective isospin Hamilto- 
nian for the doped system is of more complex type than 
the simple T=| ITF model in the undoped case due to 
the existence of biquadratic terms. Contrary to the ITF 
model there is no exact solution. However it was ar- 
gued that a mf treatment of the extended model for the 
doped ladder still gives a qualitatively correct description 
of charge order and its destruction by doping. 

Indeed it was found that the CO ground state is rapidly 
destroyed by hole doping and the size of the slope of 
the order parameter reduction -dx/d<5 increases strongly 
when one approaches the quantum critical point (A < 
1) which corresponds qualitatively to the observed be- 
haviour in NaV2 05 . A phase diagram in the doping(<5) 
- control parameter (A) plane was derived showing the 
evolution of the charge order parameter. 

Naturally the present model should be investigated be- 
yond the present mf approximation. As shown in the case 
of the ITF model for pure NaA^Os 14 the exact diagonali- 
sation Lanczos approach of finite ladders or clusters with 



the Trellis lattice structure is a powerful method to gain 
more insight into the appearance of charge order. It can 
also be used for the present extended T=l isospin model 
of hole- doped NaX^Os . 
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